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1 INTRODUCTION 



The generation and amplification of magnetic field in 
astrophysical systems are ubiquitous phenomena. The 
origin and exact mechanism of growth of magnetic fields 
in stars or galaxies have been long standing problems. It 
has been thought that helical kinetic motion or turbulence 
amplifies the magnetic field(B field). However, the helical 
component does not seem to be an absolute necessity for 
the amplification of large scale magnetic field. In astro- 
physical dynamos, for instance, the kinetic energy of some 
celestial objects like supernovae or galaxy clusters has low 
or practically zero level of helical component. The evolution 
of B fields in these objects is thought to be dominated 
by small scale dynamo(S'5'-D): the amplificatio n of fields 
below t he large scale eddy without helicitvlKazan tsev 



ABSTRACT 

To investigate the effect of energy and lielicity on tlie grovifth of magnetic field, helical 
kinetic forcing was applied to the magnetohydrodynaniic(MHD) system that had a 
specific distribution of energy and helicity as initial conditions. Simulation results 
show the saturation of a system is not influenced by the initial conditions, but the 
growth rate of large scale magnetic field is proportionally dependent on the initial large 
scale magnetic energy and helicity. It is already known that the helical component of 
small scale magnetic field(i.e., current helicity (j ■ b)) quenches the growth of large 
scale magnetic field. However, (j • b) can also boost the growth of large scale magnetic 
field by changing its sign and magnitude. In addition, simulation shows the nonhelical 
magnetic field can suppress the velocity field through Lorentz force. Comparison of the 
profiles of evolving magnetic and kinetic energy indicates that kinetic energy migrates 
backward when the external energy flows into the three dimensional MHD system, 
which means the velocity field may play a preceding role in the very early MHD 
dynamo stage. 

approximation(_ED(5Ai'M, iPouguet et al.l (|l976l )) can be 
used to explain the influences of ICs on the profile of grow- 
ing B field qualitatively. However, it is partial and incom- 
plete. Development and verification of the theoretical results 
with more detailed simulation data are necessary. Nonethe- 
less, related simulation results still provide us many detailed 
phenomena that are helpful to understanding the MHD tur- 
bulence. In this paper the effects of initial magnetic energy 
and helicity on the large scale dynamo were investigated 
using simulation data and analytic methods. 



1 196811 , iKulsrud fc AndersonI (Il992l). Menegu z zi et al. 



I 198II). iHaueen et al.l ||2003 ). ISchekochihin et al.l ll20oi) ). 2 PROBLEM TO BE SOLVED AND METHODS 

The main aim of this paper is to find out the effect 
of initial conditions(/Cs) on the growth and saturation 
of magnetic helicity(//M = 1/2(A • B), B = V x A) 
and magnetic energy (_Ejv/). For this, the combinations of 
three simulations were carried out: Non Helical Magnetic 
Forcing(iV_H'M_F), Helical Magnetic Forcmg{H MF), and 
Helical Kinetic Forcmg{H K F) . To explain simulation 
results, the equations deriv ed from EDQNM a nd tw o scale 
mean field dynamo theorv ljField fc BlackmanI (|2002r )) were 
used. 



iMininni et al.l ( 20051 )). So, it is important to understand the 
detailed mechanism of dynamo in MHD equations whether 
or not the driving force is helical. 

As of yet some problems in the MHD dynamo pro- 
cess are not completely understood: the role of he- 
lical or nonhelical kinetic(magnetic) field, the effects 
of initial conditions(7Cs) such as kinetic (magnetic) en- 
ergy and helicity. There wer e trials to see the effect s 
of ICs on t he dy namo (i Haugen fc Brandenburg (|200J), 
iMaron et al.l (|2004l )). However, the trials are not yet 
enough; moreover, there are few analytic studies to ex- 
plain the effects of initial conditions. Some statistical 
methods like Eddy Damped Quasi Normal Markovian 



For t he simulation code, high order finite difference Pencil 
Code l|Brandenburg l|200ll )) and the message passing inter- 
face(MPI) were used. The equations solved for HKF in the 
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code are, 

Dp 
Dt 

'm 

dA 
dt 



E„ and IHi,l (t, = i/ = 0.006) 



-pV- V 



(1) 



^c^Vlnp + i^ + u{V\ + ivV • v) + f (2) 
pi 



— - = vxB-r?J 



(3) 



p: density; v: velocity; B: magnetic field; A: vector poten- 
tial; J: current density; D/Dt{— d/dt + v • V): advective 
derivative; r}: magnetic diffusivity; v{=^/ p, /i: viscosity, p: 
density): kinematic viscosity; c^: sound speed; f: forcing 
function (helical or nonhelical) . The unit used in the code is 
'cgs'. Velocities are expressed in units of Cs, and magnetic 
fields in units of (popo)^ Cs (i.e., B — ^po fiov, po is 
magnetic permeability and po is the initial density). Note 
that Po ~ p in the weakly compressible simulations. These 
constants Cs, po, and po are set to be '1'. In the simulations 
?7(=c^/47rcr, a: conductivity) and v are 0.006. 

In case of the magnetically driven simulation(magnetic 
forcing, MF), forcing function f is located in the magnetic 
induction equation(9A/9t = vxB — TjJ + f) instead of 
the momentum equation. As Ohm's law(?7J = E + v x B) 
implies, f symbolizes sort of the e xternal electromagneti c 
force tha t drives the magnetic eddv ljEinaudi fc Vellil ([1999|), 
iPark fc B lackman (2012b)). 

We employ a cube like periodic box of spatial volume 
(27r)^ with mesh size of 256^ for runs. The forcing function 
f( http://pencil-code.nordita.org) used in the simulations is 
either fully helical(in fourier space, V x f = fc f , fc: wave 
number) or non-helical(V x f 7^ fcf). ({x,t) is represented 
by N {k{t) exp[i\i.{{t) ■ x -|- i(j>{t)]{N: normalization factor, 
kf(t): forcing wave number). And to prevent the shock phe- 
nomenon, forcing magnitude fk is 0.07 for KF and 0.01 for 
Af_F(note that V x f = fc/f for the helical forcing). This 
makes mach number(=i;/cs) less than 0.3. 



3 SIMULATION RESULT 1 

Fig. 1 includes the early time profiles of large scale | Hm \ (solid 



line) and _Bm (dotted line). In Fig 1(a) top line group in- 
cludes \Hm\ and Em for the case of NHMF -5> HKF: after 
Non Helical Magnetic Forcing(as a precursor simulation, 
kf=30 t<10.6) Helical Kinetic Forcing(fc/=5, t>10.6) was 
done over this preliminary simulation. The middle lines 
are \Hm\ and Em for HMF ->■ HKF: Helical Magnetic 
Forcing(fc=30, f<13.0) HKF{kf=5, t>13.0). The lowest 
lines are \H]\i\ and Em for HKF system at kf—5 as a 
reference simulation. The spectra show NHMF — > HKF is 
the most efficient in the growth of large scale Hm and EAt 
in the early time regime. HKF of which energy transfer 
chiefly depends on a effect appears to be the least efficient 
in energy transfer. In case of HMF — )■ HKF, the efficiency 
is between 'NHMF -i> HKF' and 'HKF\ And during 
HMF, the features of NHMF and HKF due to the hehcal 
and nonhelical field are observed. 

The spectrum of magnetic energy Em is always positive, 
but the sign of Hm is influenced by the external driving 
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(c) |-ff]\/| and Em (Linear scale) 

Figure 1. Seed energy and helicity in each case are the same. 
Precursor simulation {N)HMF changes the given seed field into 
the specific energy distribution, which is used as new initial con- 
ditions for the consecutive main simulation HKF. (a) Prelim- 
inary simulation NHMF{\fk\ = 0.01 at kf = 30) finishes at 
t = 10.6. During this time regime, Hm is negative. In contrast, 
HMF{\fk\ = 0.01 at kf = 30 for t < 13.0) generates positive 
Hm- HKF{\f}i\ = 0.07 at kf = 5) follows these preliminary sim- 
ulations. And HKF without a precursor simulation was done 
separately as a reference, (b) Except NHMF, the magnetic fields 
in the other cases are indistinguishably small, (c) The left line 
group shows the influence of ICs. The difference in the onset po- 
sitions is mainly decided by large scale Emo and Hmo generated 
by the precursor simulations. And right line group includes the 
shifted Ej^i and Hm of each case for the comparison. 
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(c) Ekin and Emag 

Figure 2. NHMF{fhm = 0, /o = 0.01, fe/ = 30) finishes at 
t = 10.6, and then HKF{fhk = 1, /o = 0.07, kf = 5) begins. 
Initially, only tiny Em is given(i?fci„ is zero). But Eki„ grows 
quickly, catches up with Em till t ~ 0.2, and outweighs it. (a) 
i?fci„ which is transferred from magnetic eddy through Lorcntz 
force migrates backward and forward, (b) The diffusion of energy 
among magnetic eddies without a effect in NHMF is tiny. Except 
the forced eddy, the energy in magnetic eddies seems to be mostly 
from kinetic eddies. After the precursor simulation, the peak of 
i?M (nonhclical) at A; = 30 disappears within a few time steps, (c) 
Comparison of Efcin ^i^d Em- 



(c) Ekin and Emag 

Figure 3. HMF{fhm = 1, /o = 0.01, kf = 30) finishes at 
t = 13.0, and then HKF{fhk = 1, /o = 0.07, kf = 5) begins, 
(a) Ekin of HMF is smaller than that of NHMF. (b) Emag of 
HMF is also smaller than that of NHMF. The second small peak 
around fc = 9, 10 is the inversely cascaded energy due to a effect. 
This peak moves backward to be merged into the new forcing 
pcak(fc = 5) when HKF begins. The peak of Em at kf = 30 also 
disappears within a few time steps. 



function '/' and forcing method. For example, in case 
of HMF, Hm and '/' have the same sign(middle lines, 
t <13.0), but in case of HKF, Hm has opposite sign of '/'. 
The cusp in this group(t ~ 13) is the rapid change of Hm 
from positive(ffMF, fhm = {ka-h) / {h'^)=+l, full helical) 
to negative(i:fis:F, /fefe = (v-cj)/(fc v^)=+l, fuU helical). The 
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Figure 4. (a), (b) are the preliminary simulation before HKF. 
(c) This plot is the same as that of (a), but the forced eddy is 
kf =5, closer to the large scale. It shows basic profile of i?fci„ 
does not so much depend on the position of forced eddy. Also 
linear and uniform kinetic energy distribution implies the energy 
transfer is more local and contiguous rather than nonlocal. 



reasons of reversed or equal sign of large scale magnetic 
helicity in HKF and H M F were partially e xplaine d in 
IPark fc BlackmanI l|2012ah . IPark fc BlackmanI (J2012lJ ). In 
contrast, the magnitude of \Hm\ tends to be continuous, 
which implies the relation between Hm and Em. This 
will be dicussed again. On the contrary, the direction and 
magnitude of Hm in NHMF are irregular because of the 
fluctuating '/'. During NHMF{t<,10.Q) the actual sign of 



Fig 1(b) includes the linearly scaled plots in Fig l(a] 
simulations started with the same seed field. 



All 



However, 

as the plots show, the system driven by NHMF has the 
largest and fastest growing Em and Hm in the early time 
regime. In addition, during HKF after the preliminary 
NHMF, Em and Hm still grow fastest, and maintain the 
largest values until they get close to the saturation. This 
shows the effects of ICs and forcing method clearly. 



Fig 1(c) is to compare directly the onset and saturation 
of Em and |i/M|. The group of three lines in the left 
part has the original data plots, and the lines in the right 
group are their shifted plots for comparison. The onset 
point of NHMF ->• HKF is the earliest, followed by 
HMF -^ HKF, and then HKF. The order of onset 
position is closely related to the different amount of mag- 
netic energy and helicity generated during the preliminary 
simulation. However, as Table. 1 shows, Hm,l{L: large 
scale), Em^l, Ekin,L, and Ekin,s{s: small scale) in each 
case are of similar values at their onset positions in spite 
of the different ICs. They are sort of critical values for 
the onset, and the time to reach this critical point is 
inversely proportional to the magnitude of ICs. Since the 
forcing method after the preliminary simulation is the same 
(HKF), ICs are determinants of the onset time. However, 
the saturation of turbulence dynamo is independent of ICs. 
Rather the saturation is decided by the external forcing '/' 
and intrinsic properties of system like viscosity i/(~_Re~^, 
Reynolds number Re — VrmsL/u. L: characteristic linear 
dimension) or magnetic diffusivity r){^ReM~^ , Magnetic 
Reynolds number Rbm = VrmsL/r]). All three simulations 
have the same saturated i?eM~30. Those features of critical 
values and saturation imply there is no long lasting memory 
effect in turbulence. This validates Markovianization in 
MHD equations for closure. 



Fig |2(a)[^2(c)| are the spectra of Ekin and Em of 
NHMF -5> HKF. During the preliminary NHMF 
at kf — 30, most Em is localized at the forced eddy 
(Fig|2(b)|. But Ekin which is larger than Em in most 



range spreads over from large to small scale(Fig 2(a) 2(c) I 



This relatively linear profile of Ekm spectrum indicates the 
pressure that makes the system homogeneous transporting 
the energy forward and backward is dominant in the 
very early time regime. On the other hand, these figures, 
especially Fig |2(c)[ imply the relation between Ei^in and 
Em. Initially only tiny seed Em was given to the system. 
However, once NHMF started, E^in caught up with Em 
by t ~ 0.2 and outweighed it. Soon, larger Ekin gets 
to induce Em which is the source of Hm- That Ekin is 
one of the sources of E m is coincident with t he result of 
EDQNM approximation (|Pouauet et al.l l| 19761 )'). On the 
other hand, the backward transfer of kinetic energy seems 
to contradict the accepted theorem that inverse cascade 
in three dimensional magnetohydrodynamic turbulence is 
not possible. However, when the energy or vorticity is not 
conserved, Ekin can be inversely cascaded. We will come 
back to this problem later. 
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Fig |3(a)f - j3(c)] are Em and Eki„ for HMF -^ HKF. 
HMF shows two kinds of energy transports: nonlocal 
transport of Em and local transport of Ekin- The for- 
mer is caused by helical field(a effect), and the latter is 
caused by the pressure. In Fig 3(a) backward migration 
of Ekin indicates the role of pressure. And the increase 
of Em in large scale implies the direct energy transfer 



from kinetic eddy(Fig 3(b) 3(c) I. Furthermore, the helical 
driving source generates current helicity (j ■ b)(=(fc^a • b)) 
and kinetic helicity (v • uj), which forms a effect in the 
system. This a effect generates the secondary peak of 
Hm around k ~ 10(Fig3(b)|. It keeps moving backward 
and merges into the main peak when HKF begins at kf — 5. 



The difference between Fig 4(a) 'NHMF) and 
Fig |4(b)[ iJ"A/_F) is just magnetic helicity ratio(//im — 0, 
fh,n = 1). Fig |4(a)| and Fig^^NHMF) have the same 
helicity ratio fhm = but different forced eddies(fc/ — 30, 
kf = 5). These clearly show kinetic energy migration and 
the basic profile of field evolution do not depend on the 
position of forced eddy so much. 



Fig 5(a) ~- 5(c) include the profiles of dHAf/dt, dEjii/dt, 
Sm(x0.005), ii'A/(x 0.005), and ((v-a;)-(j-b))/2(x0.001). 
These plots are helpful to infer their relative roles and 
relations in large scale dynamo(Eq.(|4]), (O). The profile of 
((v-a;)-(j-b) ~— q) shows the effects of '/' and ICs clearly. 
With the smallest ICs in HKF{Ta\Ae.l), its duration 
time(0 < t <~ 200, Fig 5(c) I of constant a coefficient is 
longer than that of other cases. In contrast the simulation 
of 'NHMF -5> HKF' has the shortest duration time of 
constant a coefficient(~ 20 < t <~ 110, Fig 5(a) I. As an 
another feature, dEM /dt is not always larger than d\HM\/dt. 



Fig 6(a) shows the profiles of d\HM\/dt in the early time 
regime. And d\HM\/dt in Fig 6(b) are shifted plots for the 
comparison. All d\HM\/dt converge to zero; but, the pro- 
file of HKF follows different paths. d\HM\/dt of HKF is 
smaller than that of the other cases until it reaches the 
onset position. HKF starts with the smallest ICs, but all 
quantities except Em,s become the same as those of other 
cases by the onset position. When the field is about to arise, 
Em,s of HKF{^ 4.3 X 10~^) is smaller than that of other 
cases(~ 7 x 10~^). In theory, this term is discarded because 
of the se emingly little i nfluen c e on the evolution of Em 
or \Hm\ (iPouauet et all (|l976l '). iField fc BlackmanI (|2002D . 
iBlackman fc Fieldl ( 20021 )'). However, Em,s is closely related 
to the conservation of magnetic helicity in the system and 
constraining velocity fleld. We will discuss about this again. 



4 THEORETICAL MODEL 

There is no theoretical method that can completely 
explain the influence of ICs like Ekin, Em, or Hm on 
MHD dy namo yet. How e ver, some approximation like 
g£>QJVM (|Pouauet etHI ^M)), though limited, can 
be used. The representations of Hm and -Ea/ of this 
method are quite similar to t hose o f two scale mean 
field m ethod (JBlackman fc Fieldl ^2004 ). iField fc BlackmanI 
(|2002l ')'). The equations are composed of Alfven effect 
term by the larger eddies, a effect term by the smaller 



eddies, and dissipation term. These approximate equations 
assume the field is composed of helical and nonhelical 
part. If helical component in the field is zero or ignor- 
ably sina.\l{N HAIF or NHKF), these equations are not 
valid. The system is divided into large(fc — 1) and small 
scale(fc = 2 ~ kmax), and this small scale can be subdivided 
into the forcing(fc = 2 ~ 6) and smaller scale(fc = 7 ~ kmax)- 



Alfven effect 



dlh^ 
dt 



dFhi 
dt 



[ {T/k){Hv - k^Hu) + rEv)+ c^Em - 2vvk''Hj, 



a Em — 2uvk Hm, 



Alfven ef feet 



(4) 



■ kT{Ev - Em) -^VHv)+ a^k^HM - 2i/„fc^SM 



^ ak^HM - iVvk'EM {E^n = Ev) 



(5) 



The coefficients are. 



[ / eky,{t){Hv{q) - q^HM{q))dq] , (a < 1) 



fc/a 



9kpq{t)Ev{q)dq, 6kpq{t) 



k/a 



1 - exp{- flkpqt) 

f-l'kpq 



Atfe = C,[l q^{Ev{q) + EM{q))dqf''' 
Jo 

-f-(l/^/3)fc[2 / EM{q)dqY^'' + {v + ri)k'', 
Jo 



r = -fc / Okpqit)EM{q)dq, 



9kpq{t)q^HM{q)dq. 



(6) 



_ff„ is kinetic helicity(=l/2(v-a;)), Hm is magnetic 
helicity(=l/2{A ■ B), please note the coefficient), and 
Vv is kinetic eddy diffusivity. a^ that transfers Hm 
and Em to larger scale is composed of the residual 
helicity (g^_ffA/(g) — Hv(q)) and triad relaxation time 9kpq- 
Ohpq is the function of eddy damping rate(see appendix) 
l^kpq{= A'fc + A^p + M?)] and connects smaller scale eddies 
and larger scale eddies. 

The infiuence of Alfven terms(fc — 0) on the large scale(fc = 
1) is physically meaningless. Ignoring Alfven terms(fc = 0), 
we find those coupled equations have two normal coordi- 
nates: 'Em + Hm' and 'Em — Hm'- 

d{EM + Hm) 



dt 



djEM - Hk 
dt 



= [a''~2Pv)[EM + HM), 



-{a'^ + 2vv){EM-HM). 



(7) 
(8) 



Assuming /q (a — 2vv)dt = (a — 2v)t, the solutii 



HM(t) = -[HMo{e 



(a-2v)t , ~{a + 2u)t 



, rn / (a — 2iy)t — (cK-\-2iy)t\ 



EM(t) = -^[EMo[e^ ' +e ^^ ') 
+HMo[e^ ' -e "- ^ ' )\. 



(9) 



(10) 



6 Kiwan Park^ E-mail: pkiwan@pas.rochester. edu^ Department of Physics and Astronomy, University of Rochesi 



These solutions show how Hm and Em are generated. For 
example, in case of HM(t) both Hmo and Emo are sources 
of HM{t), but Emo produces HKiif) like an auxiliary 
source(~Sinh). In the early time regime the effect of Emo 
on HM(t) is tiny, but finally becomes on a level with 
Hmo- HM,sat converges to -BM,aat(Fig|f (c)l) as i-^oo. This 



S.OxlO"' 



Nonhelical MF(tS10.6, k, = 30) ^ Helical KF(k[=5) 



solution shows that as long as a is larger than dissipation 
2zJ, large scale mag netic field eventually becomes fully 
helica l by a eflFect. (|Stribhng fc Matthaeuj (|l99j), iBiskampI 

liooi)) 

On the other hand, one of our interests is how long the 
effects of initial values last in turbulence. In iMaron et al.l 
(|200J) , the influence of imposed large scale magnetic energy 
on the system was tested(nonhelical kinetic forcing). The 
strong magnetic held in the large scale was expected to 
suppress the formation of small scale flelds. However, the 
effect of imposed magnetic energy disappeared soon, and 
the system eventually followed the external forcing source. 
Like an oscillator driven by an external driving source, the 
effect of ICs exist only in the early transient mode. In 
turbulence smaller eddy loses the information faster than 
large eddy does. 

To see how long the effect of ICs lasts, the formal solution 
of Eq.Q may be useful: 




(a) 

Helical MF(tS13, k,=30) ^ Helical KF(k,=30) 



Hi 



r(t„)=e ' "' ' [la'' ""- ' a (T)EM{T)dT 
Jo 




-Hm{0)]. (11) 



Using the trapezoidal method for the integration part with 
the assumption of J " 2!^i,(r)dr = 2V{tn) and t„ = nAt, we 
find the approximate solution: 



(b) 

Helical KF(taO, k[=5) 



Hm {tn 



2V(0)-2V{t„)R 



{0)EM{0)At- 



-2V(t„) 



HAf(O)] 



r 2V(ti)-2V(t„) R/, ■, rp /, \ , 2V(t2)-2V(t„) Ru ^ c^ u \ , 

+ [e '■^' "-"'a (ti)EM(ti) + e '^' '"'a (t2)EM{t2) + 

I „2V(t„_i)-2V(t„)^fl/, X ri ,-+ ^^ A + 

These show all previous results affect the current magnetic 
helicity in principle. However the influence decreases expo- 
nentially, which is coincident with the simulation results. 
The decaying speed depends on the several factors: energy, 
helicity, v, and 77. Of course the actual Vv varies with time. 
But, since i'„(~ V) changes rather smoothly and satu- 
rates to a constant, this inference is qualitatively reasonable. 



5 SIMULATION RESULT 2 

Table. 1 provides information on the energy distributions 
of the evolving variables at each stage. When Em or J/a/ 
is about to rise, most Ekin is located in the forcing scale 
regime(A: = 2 ~ 6) regardless of its initial distribution. 
The ratio of smaller scale Eki„{k = 7 ~ kmax) to Ekin of 
the whole small scale regime(A; = 2 ~ kmax) is about 1%. 
After the onset, as large scale Em or Hm grows, kinetic 
energy migrates towards the smaller scale. At this time 
the saturated ratio elevates up to 8 ~ 10%. When large 
scale eddy needs more Hm or _BAf (onset position), more 
(helical) kinetic energy is located in the forcing scale. And 
if the inverse cascade of Hm is less required, Ekin in the 




dlH„l/dt 

dFVdt 

IH„l(xO,005) 

hjxo.005) 

(vc.-i-b)/2(x0.001) 



200 300 

time 



400 500 



(c) 

Figure 5. During HKF, {(v-lj)-(j- b))/2 in each case drops at 
different time position. It depends on the energy and helicity(7Cs) 
from the preliminary simulation. Also due to the different eddy 
turnover time between large and small scale, there is a phase dif- 
ference in the profile of growth rate, Em{Hm), and a related 
term. For the growth ratio, usually logarithmic growth ratio is 
used: dlog\HM\/<it = —oi^Em /\Hm\ — 2k'^Uy (k = 1), but linear 
growth rate was used for the mathematical convenience and visi- 
bility. All quantities but Em a^rid Hm are the averages of 50~ 100 
nearby values. 
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— dlH„l/dt NHMK(k,= 30. IS 10.6)^HKF(k,= 5) 
- dlH„l7dt HMF (k,= 30, IS 13,0)^HKF(k,= D)l 
- dlH„l/dt HKF(k, = 5, tgO) 



(a) 
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1.0x10"*- 



dlH„l/dt (k=l) 
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HMF^HKF 

HKF I 




200 300 

time 



(b) 



Figure 6. (a) Growth ratio proportionally depends on the ICs. The area between the line and time axis is Hm- (b) The profile of 
evolving growth ratio of HKF is slightly different from that of others. It seems to be caused by the turbulent effect. Each growth ratio 
is the average of 50 nearby points. 
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Figure 7. The direction of magnetic helicity is decided according to the conservation of total magnetic helicity in the system. The 
minimum {t~250) and turning point(t~280) of magnetic helicity in (b) can be compared with the change in growth of large scale Ht,[ 
in Fig|5(^ 



forcing scale decreases and moves toward the smaller scale 
which has less helical effect but more dissipative effect. 
Ekin,3, more exactly (v • uj) plays the role of balancing 
the growth of large scale magnetic field. However, the 
evolution of Ekin,L shows rather an irregular feature. In the 
early time regime Ekin,L leads the growth of Em^l which 
generates HM{¥ig.l, 2, and 3). According to EDQNM 
approximation, the role of Ekin,L with Em,l is related to 
the self distortion effect (the eddy damping rate ^t, Eq.(|6])). 
Smaller magnitude of Ekin,L decreases fik, which increases 
a effect and dissipation at the same time. However, more 
detailed simulation is necessary to check these theoretical 
inference. 



Analytic equation like EQDNM or mean field dynamo 
theory does not explicitly explain the role of magnetic 
energy in the small scale. But simulation results provide 
some clues to the influence of Em^s on the dynamo. The 
ratio of smaller scale Em to that of the whole small 
scale is consistently regular, i.e., from onset: ~ 30% to 



saturation: ~ 20%. The distribution of Em, 3, more exactly 
(j • b), is related with the inverse cascade of Hm {01 Em) 
and balancing the growth rate of large scale magnetic 
fleld. {j ■ b) in a coefficient does not always quench the 
large scale magnetic field. As Fig 7(a) 7(b) show, when 
the necessity of inverse cascade of magnetic energy is 
large, forcing scale Hm is negative so that a effect is en- 
hanced(the kinetic helicity in a coefficient keeps positive). 
As the large scale field saturates, the sign of forcing scale 
magnetic helicity grows to be positive, i.e., lowering a effect. 



With more detailed plots we can investigate the dynamic 
properties of small scale regime with respect to the large 
scale B field growth. 

shows the evolution of small scale 
2 EM{k)) and nonhelical Em,s{~ 
Em,s — kHM,s, thin line). To remove the preliminary 
simulation, NHMF — > HKF{= 1st simulation from now 
on) was shifted by -10.6 time unit using t' = t — 10.6 and 
HMF -^ HKF(= 2nd simulation) was shifted by -13.0 



Fig 8(a) 



77 / sr-^k — k 
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2.7 X 10"^ 



1.2 X 10"" 



6.5 X 10"" 
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-5.7 X 10" 



-4.9 X 10" 
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TT" 



'^M,L 



7.7 X lO"'' 

(2.5 X 10"^) 

(32%) 



5.7 X 10-^ 

(1.2 X 10"^) 

(21%) 



5.5 X 10"^ 

(5.5 X 10"2) 

(100%) 



2.7 X 10" 



5.8 X 10"^ 

(1.1 X 10"^) 

(19%) 



2.4 X 10" 



6.6 X 10" 



2.7 X 10" 



4.3 X 10"-'= 

(1.5 X 10"^) 

(~35%) 



Em,s 



5.6 X 10"^ 

(5.6 X 10-2) 

(100%) 



2.8 X lO""" 

6.5 X 10-^ 

(6.3 X 10"^) 

(97%) 



1.7 X 10""^ 

2.1 X 10-2 

(1.8 X 10-*) 

(1%) 



1.1 X 10-"^ 
9.7 X 10-^ 



1.6 X lO"**" 



7.5 X 10- 
(2.5 X 10" 
(33%) 



1.7 X 10" 



8.6 X 10" 



1.7 X lO-'* 

(1.7 X 10-**) 

(100%) 



3.3 X 10-^"~ 

1.3 X 10-^ 

(2.7 X 10-*) 

(- 0%) 



1.9 X 10~^ 

2.1 X 10-2 

(2.4 X 10-"') 

(1%) 



6.1 X 10 
(1.3 X 10 

(21%) 



4.3 X lO--'' 

1.1 X 10-2 

(8.5 X 10-"') 

(8%) 



Ek.l 



7.7 X lO-' 

(6.8 X 10-^) 

(88%) 



2.3 X 10-2 

(2.7 X 10-") 

(1%) 



9.5 X 10-^ 

(9.2 X lO-"*) 

(10%) 



Ek,s 



(8.6 



X 10-*) 
(9%) 



Table 1. Large scale: A: = 1, forcing scale: fc = 2 ~ 6, smaller scale: k 
The percentage is the ratio of smaller scale to small scale:~ J^fe™?"" 



/Ej" 



^(quantity in parentheses). -EAf,a(*) = T^.i. 



,{k,t). 



time unit using t' = t — 13.0. But HKF{= 3rd simulation) 
was not shifted. The difTerence between the thick and thin 
line is the helical B-field, which is generated counteractively 
as the large scale helical magnetic field grows. This helical 
magnetic field in the small scale can induce the growth of 
large scale B-field({j ■ &) > 0) or suppress it{{j ■ h) < 0). 



Fig 8(b) shows the profiles of small scale Em in the early 
time regime. Em of the 1st simulation is relatively larger 
than that of other cases, and Em of the 3rd simulation is 
the least. The initial small scale Emo,b{5.Q x 10-2 of the 
1st simulation, 5.5 x 10-2 of the 2nd simulation) due to 
the preliminary simulation drops till i ~ 10 and begins to 
grow again. The different minimum value and evolution of 
each field profile imply some important clues to the relation 
between large scale magnetic field and small scale magnetic 
field. In addition, the origin of helical magnetic field can be 



inferred from Fig 8(a) 8(b) 



The magnetic helicity is the topological linking number of 
magnetic fields. But statistically it can be considered as 
the correlation betwe en different components of magnetic 
field (|Yoshizawal l|201ll )). 



{B^B.i-k)) = P^Ak)^^ + ~e^,iHM{k) 



(12) 



Fij \k) — dij 



"fc2" 



Since magnetic helicity (a ■ b) cannot be larger than 2Em /k. 
Small scale magnetic energy has a lower bound proportional 
to the small scale magnetic helicity. In addition, the growth 
of small scale magnetic helicity depends on that of large 
scale magnetic helicity in terms of the conservation of 
magnetic helicity in the system. All of these explain the 
reasons of quick drop of Emo and different evolution of 
Em in the small scale. The role of nonhelical Em becomes 
clear with the comparison of i5fci„. We will discuss about 
this again. The initial B-field plays the role of seed field 
in MHD dynamo, and at the same time the correlation 
between its different components constrains the growth of 
large scale magnetic field dynamically changing the sign 
and magnitude. 



Fig 8(c) 8(d) show Ekin spectra of the 1st and 2nd 
simulation are very similar. Fig |8(b)[ |8(c)[ and |8(d)| imply 
the profile of kinetic energy does not depend on the 



small scale Em much as long as Em is not too much 
different. 



Fig 



and 8(f) also show large scale Ek 



is independent of Em when the magnetic energy is not 



significantly different. In addition. Fig 8(d) 8(e) clearly 



show that Ekin drops when large scale Em begins to rise, 
i.e., onset position. At the onset point of the 1st simulation, 
t ~ 70, Ekin of this simulation begins to drop. And 
around t ~ 120, onset position of the 2nd simulation, Ekin 
of this simulation also begins to drop. They meet again 
each other when large Em of each simulation gets saturated. 

In Navier Stokes equation(Eq.2), Lorentz force(J x B) 
can be decomposed into magnetic tension(B-VB) and 
pressure(— VB2/2). The force parallel to B-field from 
magnetic tension and pressure is canceled out. Only 
the force perpendicular to the magnetic field li ne like 
(B'^/R c - V(B2/2))n, Re. radius of curvature, IPriestI 
(|2003i )) exists, as the definition of Lorentz force implies. 
When the growth of Em accelerates near the onset point, 
the compressive force(— VB2/2) normal to B-field grows 
so that the net effect of Lorentz force becomes negative. 
This presses the plasma and causes the geometrical changes 
of magnetic fields. The kinetic motion of plasma slows down. 



Figure 8(c) 8(f) also show there is a time regime(t < 7 ~ 10) 



where the profile of evolving Ekin is independent of the 
initial values. This occurs when Lorentz force is still weak, 
and looks like the corresponding concept of the kinematic 
regime. 



Fig 9(a) ~- 9(f) include kinetic and current helicity. Kinetic 
helicity profiles shown in Fig 9(a) are very similar to 
those of kinetic energy. In addition, Fig |9(b)| shows the 
clear migration of kinetic helicity toward the smaller scale 



regime. Fig 9(c) 9(d) are the current helicity in small 



scale regime. Especially the profile of evolving current 
helicity in Fig |9(d)| suggests that the evolution of Em in 
small scale be determined by the growth of large scale 
magnetic energy (helicity). As mentioned, small scale Emo 
of the 1st and 2nd simulation are almost the same. We 
know the different large scale Emo(ox Hmo) causes the 
different growth rates(9B/(?i ~ (v • ti;)B). The fast growth 
of negative Hm in the large scale requires the fast growth of 
positive Hm in small scale to conserve Hm in the system. 
The evolution of small scale Em or Hm is highly influenced 
by the initial large scale Em or Hm- The quick change of 
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Figure 8. Preliminary simulation effect is removed by shifting the time unit(lst simulation: -10.6, 2nd simulation: -13.0). (a) Emo.s- 
5.6x 10~'^(lst), 5.5x 10~^(2nd), 1.7x 10~*(3rd). (c), (f) In the very early time regime, all B^i^s evolve in the same way being independent 
of the initial values, (d) i?fei,j of smaller scale eddy branches off earlier. 



sign with the sequent fast growth of small scale (j • b) of 
the 1st simulation and the fast decay of small scale Em of 
the second simulation support this fact very well. 



Fig |9(e)[|9(f)| show how residual helicity evolves. In the very 
early time regime the residual helicity does not depend on 
the initial conditions. This phenomenon is an inevitable 
consequence for the kinetically driven MHD dynamo. 



In Fig 10(a) ^ 10(f) the onset positions of 2nd and 3rd 
simulation were shifted by i — > t — (120 — 70) and 
t ^ t — (186 — 70) to compare the behaviors of field profiles 
after the onset. Except some minor differences due to 
the turbulene, all three simulations have the same field 
profiles. This indicates the saturation of MHD dynamo is 
independent of the initial conditions. 



On the other hand in Fig 4(a) ^ 4(c) we pointed out the 
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Figure 9. Small box includes (v ■ u)) and residual helieity for 3rd simulation. Their profiles are influenced by the large scale field of the 
3rd simulation and take different routes unlike 1st and 2nd simulation, (d) The evolution of (j • b) needs to be compared with that of 
Em.s in Fig|8(b)| 



cascade of kinetic energy toward larger scale. This looks to 
contradict the established theory that the inverse cascade 
occurs with the strong rotation effect or in the ideal (quasi) 
two dimensional hydrodynamic system. This conclusion is 
based on the conservation of physical quantities l ike energy 
and enstrophy. The modified expression of ik) (JDavidsonl 



(J2004 l'l is, 
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Figure 10. The plots of 2nd and 3rd arc shifted by -50 and -116 for the comparison with 1st simulation after the onset. 



d{k) _ d n^kE^{k,t)dk ' 
dt ^ dt\ j^ E^{k,t)dk 

-2ki[irik-k^fE4k,t)dk] 
J^ E^{k,t)dk 

2fci[/o°°('='+^/)-E.(fc,i)rffc] 



+ 



j^ E^{k,t)dk 
{!^kE^{k,t) dk)j-J^ E,{k,t) dk 
ij;°E4k,t)dk)^ 



(fcc = / kE{k, t) dk/ J E{k, t) dk) 



If total energy J°° Ey{k,t)dk and enstropy 
J^ k^ Ev{k,t) dk are conserved, the first term in tiie 
right hand side determines (k). Since the usual spreading 
Ev{k,t) in turbulent flow makes J^{k — kc)^ Ev{k,t) dk 
grow, (k) decreases (inverse cascade). However, if enstro- 
phy or energy is not conserved(by the external source), 
(k) can grow or de crease according to E^ and {uj^/2). 
iBiferale et al.l (|2012l ) showed that the reverse cascade of 
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hydrodynamic energy occurs when the mirror symme- 
try is broken (heUcity). In addition, there were another 
trials to explain the forward cascade of E ki n using the 
canon ical ensemble average (|Biskamp| (|2008h . iFrisch et alj 
([1973)). If Etot, Hm, and Hci={u ■ b)) are conserved 
quantities(ideal three dimensional MHD system), the 
form of i?fcin(~ {uiUi}) calculated using the canonical 
partition function Z~^exp{—aEtot — PHm — "/He) {Z: a 
normalization factor) shows forward cascade. However, if 
Etot{= Em + Ekin) is not conserved because of the external 
forcing or some other sources, another term should be 
added to Etot in the partition function; this changes the 
averaged Eki,i into a new form that allows the backward 
cascade. 



6 CONCLUSION 

Based on simulation and theory, we have investigated the 
influence of ICs, the role of (non)helical field with the 
pressure in the energy transfer. 

The growth rate of large scale magnetic field is chiefly 
proportional to the large scale initial values. In contrast, its 
saturation depends on the external driving source and the 
intrinsic properties like -q ox v instead of ICs. 

Comparing the simulation results, we have seen how the 
helical and nonhelical magnetic field constrain MHD dy- 
namo. The helical magnetic field in the small scale has been 
thought to quench the growth of large scale magnetic field. 
As (A ■ B) increases the opposite sign of small scale (a ■ b) 
also grows. So the amplification of large scale magnetic 
field slows down and saturates. However, in the early time 
regime (a ■ b) and large scale Hm have the same sign; thus, 
(a • b) boosts the growth of large scale magnetic field. 

On the other hand, growing nonhelical magnetic field 
presses the plasma through Lorentz force(magnetic pres- 
sure) and slows down the motion, which constrains the 
magnetic fields eventually. Of course, the evolution of small 
scale fields is also influenced by the large scale ICs and the 
evolving large scale field. 

Besides, it is observed kinetic energy migrates backward 
when the external energy flows into the three dimensional 
MHD system. And the velocity field in the early time regime 
seems to play a preceding role in the MHD dynamo. 
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APPENDIX A: EDDY DAMPED QUASI 
NORMAL MARKOVIANIZATION 



Al Tavo point closure 

Navier Stokes equation for the incompressible fiuid is 

d_ 
dt 



— +uk'^jui{k) = -ikmPijik) ^ Uj(p)um(g) ~ (uu). (Al) 

p+q— fc 



This equation requires information on the second order cor- 
relation equation. Then we need to solve another differential 
equation: 

d_ 

dt 



+ i/(fc^ + k'^)^ {u(fc)u(fc')) = (uuu). 



(A2) 



We can derive the third order correlation term, which needs 
the fourth order correlation: 

^ + i^ik^ + / + q^)\ (u(fc)u(p)u(g)) = (uuuu). (A3) 

It is known that the probability distribution of turbulent 
velocity is not far from the normal distribution. Then, the 
fourth order correlation term can be decomposed into the 
combination of the sec ond order correlat i on te r ms(Quasi 
Norm al approximation, IProudman fc ReidI (Il954r ). iTatsumil 
|l953)). 



{u{k)u{p)u{q)u{r)) ~ ^(uu)(uu). 



(A4) 



A2 Eddy Damping coefficient 

However, lOgural (| 19631 ) pointed out that Quasi Normal 
appro ximatio n coul d make the energy spectrum negative. 
Later lOrszag (|l97Ct ) found that the decomposed value be- 
came too large when the fourth correlation was decomposed 
of the combination of second correlation terms. Orszag in- 
troduced eddy damping coefficient /ifcpq(~ 1/i)- 



— + u(k'^ +p^ + q^) + ^J.kpq j (u(fc)u(p)u(g)) 
Orszag suggested 

^lk ~ [fc'£;(fc)]l/^ 



^J.kpq = /ife + Mp + Mqi 



(A5) 



(A6) 



(Mfcpij used in Eq[S]is a little different from Orszag's one.) 
However, if energy drops faster than fc""^, eddy damping 
term(~ i~^) decreases with 'fc'. This means the damping 
time of a smaller eddy can be larger than that of a larger 
eddy. To solve thi s problem, another modified representation 
was suggested by iLesieur fc Schertzen (|l978l ') : 



Hk 



E{p,t)dp 



1/2 



(A7) 



This is 'Eddy Damped Quasi Normal' approximation. 
Then we have, 



dt 



[uiUj)k. 



dr I e 

Jk+p+q=0 



-[Mfepg+i'Ci: +P +9 )](t- 



'"•' y^ (uu) (uu) dp t 



(A8) 
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If time scale of 5]](uu)(uu) is much larger than [nkpq + 
u{k'^ +p'^ +q'^]^^ , markovianization makes the equation much 
simpler. 

d 



dt 



+ 2uk ]{uiUj)k,t 



I 



^kpq 



Okpq ^{uu)(uu)dpdg. (A9) 
(AlO) 



This is called ^Eddy Da mped Quasi Normal Markovian 
appro x imation' ( EDQN M. IPouauet et~al] (|l976l ). IPavidsonl 
l|200ll ). lLesieuil (|200^ )). 
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